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ABSTRACT 

Filaments, forming in the context of cosmological structure formation, are not only 
supposed to host the majority of the baryons at low redshifts in the form of the warm- 
hot intergalactic medium (WHIM), but also to supply forming galaxies at higher 
redshifts with a substantial amount of cold gas via cold steams. In order to get insight 
into the hydro- and thermodynamical characteristics of these structures, we performed 
a series of hydrodynamical simulations. Instead of analyzing extensive simulations 
of cosmological structure formation, we simulate certain well-defined structures and 
study the impact of different physi cal processes as well a s of the scale dependencies. In 
this paper, we continue our work (|Klar fc Mucketll2010() . and extend our simulations 
into three dimensions. Instead of a pancake structure, we now obtain a configuration 
consisting of well-defined sheets, filaments, and a gaseous halo. We use a set of sim- 
ulations, parametrized by the length of the initial perturbation L, to obtain detailed 
information on the state of the gas and its evolution inside the filament. For L > A 
Mpc, we obtain filaments which are fully confined by an accretion shock. Additionally, 
they exhibit an isothermal core, which temperature is balanced by radiative cooling 
and heating due to the UV background. This indicates on a multiphase structure for 
the medium temperature WHIM. We obtain scaling relations for the main quantities 
of this core. After its formation, the core becomes shielded against further infall of 
gas onto the filament, and its mass content decreases with time. In the vicinity of the 
halo, the filament's core can be attributed to the cold streams found in cosmological 
hydro-simulations. They are constricted by the outwards moving accretion shock of 
the halo. Thermal conduction can lead to a complete evaporation of the cold stream 
for L > 6 Mpc/h. This corresponds to halos more massive than Mhaio = 10 13 M Q , 
and implies that star-formation in more massive galaxies can not be supplied by cold 
streams. 
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1 INTRODUCTION 

At any time during the cosmic history, the majority of 
baryons is supposed to rest not in stars and galaxies, but 
in the intergalactic medium (IGM). During the evolution 
toward low redshifts, this reservoir undergoes a substantial 
change in properties. At redshifts of « 2 most of the IGM is 
photoionized and is at a temperature of about 10 4 K. This 
gas can be detected in absorption in the spectra of high red- 
shifted quasars. At low re dshifts, only about 20 % of th e 
IGM can be in this state IjProchaska fc Tumlinsonll2009f) . 
The rest can not be accounted for by observational means. 
The explanation for the absence of these missing baryons, 
which is currently favored in the community, is connected 
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to the formation of the luminous structures in the Universe. 
While streaming into the gravitational potential wells of the 
forming structures, the velocity of the gas eventually be- 
comes supersonic. This leads to the formation of shocks, 
which are able to convert the kinetic energy of the gas 
into thermal energy. By this shock heating, the mean-scale 
streaming motions produce shock-confined sheets and fila- 
ments which contain gas at much higher temperatures com- 
pared to the photoionized high redshift IGM. This phase of 
the IGM is called warm-hot intergalactic medium (WHIM) 

The above scenario was developed using hydrody- 
na mic simulations o f cosm ic structure formation. Results 
by ICen fc Ostrikerl (J1999T ) suggest that approximately 30 
% to 50 % of IGM at z = rest in the WHIM. 
iDave et all l| 19991 . l200ll . |201Ch confirmed these results with 
simulations of increasing resolution and examined fur- 
ther the interplay between the different phases of the 
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IGM: photoionized, condensed, warm-hot, and hot. Subse- 
quently, more and more simulations were performed, includ- 
ing additional effects like models for galactic superwinds 
l|Cen fc Ostriker] 120061) and non-equil ibrium effects in the 
chemical networks (ICen fc Fans 20061 ). Using zoom-in sim- 
ulation techniaues lDolag et al.l l|2006h were able to simulate 
a whole super-cluster filament scenario with high resolu- 
tion. Further numerical simulations (e. g., Kang et al.l 120051 ; 
iKawahara et al.l 120061 : iTornatore et al.l |2010j) with different 
numerical schemes and resolutions also consistently support 
the described scenario. Recentl y, the Over Whelmi ngly Large 
Simulations project (OWLS) (jSchave et al.l I2009T ). consist- 
ing of a set of large hydrodynamical cosmological simula- 
tions, reporte d predictions for the soft X-ray and UV metal 
line emission IIBertone et al.ll20 10a b) and the absorption by 
OV1 absorbenT ifTepper-Garcia et al.llioilh . 

Due to the high degree of ionization, the observational 
signature of the WHIM is very weak, in parti cular with 
respe c t to absorption of neutral hydrogen (|Cen et al.l 
l200ll : iRichteret all l2006al lbh . Therefore, the detection of 
highly ionized metal lines in the spectra of bright quasars 
and bl azars is muc h more promising (IHellsten et all 
19981; IPerna fc Loebl Il998l : Fang fc Canizaresl I2OO0F 



nal structure . In high resolution hydrodynamical simulatio ns 
dKeres et al.ll2005l; lOcvirk et all |200S| ; iKeres et all l2009al lbl; 



Fang fc Brvanl 1 20011). These pred ic tions l ed to a series 



of detections (Tripp et al.l [2OO0I. IpOOll: INicastro et al 
20021 ; iFang et al.l 120021: iMathur et al.l 12003J: iFuiimoto et a! 



20041; IDanforth fc Shulll 120051; IDanforth et all 120061 : 
Danforth fc ShulfeOOStlTripp et al.ll2008l ), but they are con- 



sidered to be rather tentative. A detection with sufficiently 
high s ignal-to-noise ra tio is reported by I Nicastro et al.l 
|2005al rbl). More recent. IPanforthet all (|2010l ) reported the 
detection of a number of absorbers in the low-redshift IGM 
in the spectrum toward the BL Lac object 1ES1553+113. 
In particular, they discuss a triple absorber complex at 
z ~ 0.19, seemingly of multi- phase nature, and correlated 
with a large-scale filament in the galax y distribution of 
the Sloan Digital Sky Survey (SDSS) (|Abazaiian et al.l 
2009). Additionally, several tentative detections of the 
WHIM through its metal line X-ray emission are claimed 
(Kaastra ct al. 2003; Finogu enov et al.ll2003l : INicastro et al.l 
l20ld : IZappacosta et al.ll2010l ). 

The WHIM is supposed to reside in shock-confined 
filaments and sheets of the IGM. These structures are the 
result of one-and two-dimensional collapse processes. The 
basic theory for the formation and evolution of structure 
is alrea dy well understood since the sixties of the last 
century (IDoroshkevich fc Zel'Doyichlll964l;IZel'Dovichlll970l; 
Sunyaev fc Zeldovichl 1 19721 ; IDoroshkevich fc Shandarinl 
19781 ). According to these theories, the most probable 
formation process starts first with the one-dimensional 
collapse. Only if the underlying dark matter-distribution 
enters the first caustic, multi-streaming of matter leads to 
the formation of two-dimensional filaments and eventually 
knots, which are characterized by matter collapsing in three 
dimensions. This evolution is closely followed by the gas 
distribution. 

Recent studies on galaxy formation indicate on another 
very important aspect of these structures. In the classical 
picture, galaxies obtain their baryons mainly by smooth ac- 
cretion of hot gas through the ir virial shock ijWhite fc Reesl 
1 19781 : iFall fc Efstathioul Il980h . Mergers with other halos 
bring more baryons to the galaxies and influence their inter- 



ICeverino et al.ll2010l ). however, another aspect of gas accre- 
tion is evident. For massive galaxies at redshifts of z — 2— 3 a 
significant amount of gas streams along the filament directly 
onto the galactic disk. These cold streams seem to penetrate 
the virial shock and keep a low temperature of « 10 4 K. 
This could have significant impact on galaxy formation. The 
massive supply of cold gas at an early point in galaxy forma- 
tion does not only influence the star-fo rmation in the galaxy 
|Dekel et al.ll2009l ; iBouche et al.ll2010t ). but also give an ex- 
planation for the unusual morphology of cl ump-clusters and 
chain- galaxies in Hubble Ult r a Deep Fiel d fElmegrce n et al.l 
20051; IBournaud et al.1 120071 ; iBournaud fc Elmegreenl 120091 ; 



Agertz et ijboosh. l "ls o discusse d in tn 7c^e ~f 



Lyman-a blobs at even higher redshifts JGoerdt et al.ll2010l : 
Fauchcr-G iguc re et al.l [2010J) . 

It is of principal importance to investigate the detailed 
thermodynamic state and the internal kinematics of the 
structures which, on the one hand, may hide a large frac- 
tion of cosmic baryons at low redshifts, and, on the other 
hand, may have considerable impact on galaxy formation at 
higher redshifts. Their detailed modeling, however, is highly 
demanding for computational astrophysics. The treatment 
of low-density regions in great detail is difficult. Compared 
with the numerical handling of high-density matter distri- 
butions where adaptive techniques can be applied, for low- 
density regions, higher overall particle and/or grid number 
is unavoidable for an appropriate description. In addition, 
higher resolution calls for a more detailed consideration of 
local physics, e.g., star formation, feedback, contamination 
by heavy elements, etc. Altogether the computational effort 
is immensely more complicated if considered within a cos- 
mological context. Despite of the currently available highly 
developed computational techniques, an adequate treatment 
of low-density regions is still at the limit or beyond the near 
future potential. For these reasons, in order to get informa- 
tions about filaments and sheets beyond to what simulations 
of structure formation can tell us, we choose a complemen- 
tary approach: Instead of analyzing extensive simulations of 
cosmological structure formation, we simulate certain well- 
defined structures. 

In iKlar fc Miicketl (|2010l . hereafter Paper I) we 
started our investigation with the consideration of the 
one-dimensional collapse of one perturbation with a 
given length scale L. This scenario i s common l y re- 



ferred to as cosmic pancake formation (Zel'Dovich 



1970 



1984 



Sunyaev & Zeldovich 1972; Binncy 1977; Bond et al.1 
Shapiro fc Struck-Marcelll Il985l ). We were able to simu 



late the detailed physics and thermodynamics of the one- 
dimensional collapse toward a gaseous sheet at extremely 
high resolution and investigated the impact of radiative cool- 
ing, heating by the photoionizing cosmic UV background, 
non ionization equilibrium (IE), and thermal conduction on 
the configuration. Above an initial perturbation scale of ~ 2 
Mpc the collapsing gas shocks. The density and tempera- 
ture profiles are characterized by the existence of a cold and 
dense core region which is at thermal equilibrium at a tem- 
perature of T c ~ 2 x 10 4 K. This core is built even before 
shock formation, and its properties are given by the inter- 
play between radiative cooling and the energetic input from 
the UV background. The determining variable for the prop- 
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erties of the core is the initial perturbation scale L. While 
non-IE chemistry has no significant effect on the simulation, 
thermal conduction provides a mechanism to shrink or even 
evaporate the core. 

Based on our one-dimensional simulations in Paper I, 
we extend our simulation into three dimensions. When con- 
sidering three perpendicular modes, multi-streaming occurs 
and leads to the formation of a well-defined filament and a 
halo. Though the shape of the gas distribution is extremely 
idealized, it is expected to exhibit the correct density and 
temperature profiles depending on the initial scale length, 
again. In addition, in three dimensions it is possible to study 
the large-scale gas velocities, which themselves are deter- 
mined by the gravitational potential. The temperature and 
density profiles for the three-dimensional case are expected 
to be much more complex even for our highly idealized ge- 
ometry. The idealized filament-halo scenario is also suited 
to study the particular hydrodynamical conditions in which 
accretion via cold streams occurs. In these simulations a rel- 
atively cold stream forms along the filament and propagates 
deep into the halo. The existence of a natural regulation 
mechanism for the size and existence of a cold core as found 
in the one-dimensional case might indicate similar restric- 
tions for the cold streams in three dimensions. Compared to 
cosmological hydro-simulations, our setup, in particular the 
smaller simulation box and the Jeans-length based refine- 
ment strategy, results in a significantly higher resolution in 
regions of intermediate overdensity, and the computational 
expense is drastically reduced. \j 

The paper is organized as follows: In the next section 
we will recapitulate the theoretical framework for this study, 
and describe the used code. Since most of these topics were 
already discussed in Paper I, we would like to refer the 
reader to the comprehensive description therein. In Sect.O 
we present the result of our three-dimensional simulations. 
The implications of these results on the physics of gaseous fil- 
aments will be further examined in Sect. [4] In Sect. [5] we use 
our simulations in order to investigate the different modes 
of gas accretion on galactic halos. In particular, we analyze 
the occurrence of cold streams in our simulation. Finally, we 
discuss our results and possible implications in Sect. [5] 



2 SIMULATIONS 

2.1 Theoretical Framework 

Since the theoretical framework of this study is described 
in detail in Paper I, we only review the main concepts. We 
describe the baryonic component of the Universe as a poly- 
tropic gas with an adiabatic exponent of 7 = 5/3. Its evo- 
lution is governed by the Euler-eq uations, which we use in 
their super-comoving formulation (jMartel fc Shapiro! 1 19981 ; 



1 A m odern simulation, i.e. from the OWLS project jSchave et al.l 
2009) which is especially dedicated to the description of the IGM, 
uses 512 3 particles in a box of 25 Mpc/h. This results, at an 
overdensity of p/p = 100, in an effective resolution of Ri 10 6 
particles per (Mpc/h) 3 . Our simulations use, e.g. 12 refinement 
levels in a box of 4 Mpc/h. This results, at the same overdensity, in 
an effective resolution of ss 10 9 cells per (Mpc/h) 3 , which exceeds 
the resolution realized by the OWLS simulation by 3 orders of 
magnitude. 



iDoumler fc Kneben20ld ). W e assume a ACDM c osmology 
close to the WMAP5 results JKomatsu et al.ll2009l ) with the 
parameters £7a = 0.73, f2 m = 0.27, and Ho = 71 Mpc km - 
s _1 . In our simulations we do not compute the collisionless 
dynamics of the dark matter (DM) self-consistently, but as- 
sume it to follow the same distribution as the baryons. This 
assumption is used for the determination of the gravitational 
potential only. Since the gravitational potential is relatively 
insensitive with respect to small deviations of the matter 
distribution, the dynamics of the gas is almost not affected. 

This is an obviously crude assumption for the descrip- 
tion of high density structures like galaxies where the baryon 
distribution is expected to be more concentrated than the 
DM distribution. For lower densities, the clustering of the 
DM takes place on arbitrary small scales determined by the 
initial density perturbation spectrum. In contrast, the gas 
has a finite temperature and pressure. In particular, after the 
re-ionization the heating by photo-ionization due to the UV 
background radiation leads to gas temperatures of the order 
of a few 10 4 K. The gas pres sure causes a much sm oother dis- 
tribution of the gas (see e.g. lHarford et al.ll2008l ). erasing in- 
homogeneities on scale lengths smaller than the correspond- 
ing Jeans length. However, the locally enhanced DM density, 
e.g. in small halos, decreases the effective Jeans length. This 
might lead eventually to the collapse of perturbations in the 
gas distribution on smaller scales, either initiated or fur- 
ther enhanced by thermal instability. Those small clumps of 
the gas distribution within the extended filaments cannot be 
considered within the used approximation. It is a reasonable 
assumption however, that the small clumps end up finally 
in hydrodynamical balance with the ambient filament gas. 
The thermodynamical state of the latter is defined by the 
UV background radiation and the shock-accreted gas. 

Thus, we compute the total matter density by p to t = 
p/Jb, where Jb — £lb/£lm = 0.16 denotes the baryon frac- 
tion. This is an obviously crude assumption for the descrip- 
tion of high density structures like galaxies, but it is suffi- 
ciently accurate at low to intermediate densities. 

We follow a chemical network consisting of the primor- 
dial species H I, H 11, He I, He 11, He III, and electrons in ioniza- 
tion equilibrium. For the computation of the chemical source 
term, as well as radiative cooling and heating by the cosmic 
UV background, we include contributions from photoion- 
ization, collisional ionization and recombination, dielectric 
recombination of Hen, collisional excitation of Hi and Hell, 
and bremsstrahlung (three-body process es are neglec ted). 
The corres p ondin g rates are taken from iBlacld l|198ll ) and 
iKatz et al.l |l996|) and are summarized in Table A.l of Pa- 
per I. We employ primordial mass fractions of hydrogen 
Xh = 0.76 and helium xh c = 0.24. 

2.2 Code 

The simulations of this paper are carried out using the 
publicly availa ble simulation code RAMSES, introduced in 
iTevssierl J2002). RAMSES is one of the most popular codes 
in numerical cosmology. It combines an Adaptive Mesh Re- 
finement structure with an N-Body solver for the DM and a 
hydrodynamic solver for the baryons. The main application 
of RAMSES lies in th e simulation of larg e -scale stuctures and 
galaxy formation jlOcvirk et al.l 12 008: Agcrtz et al. 2009, 
l201ll : iGoerdt et ail 120101 ; iTevssier et all l2010l : lHahn et al.l 
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2010). It is, however, also applied in simu lations of star- 
form a tion and inter-stellar medium physics ( Froman g et al 



20061 : iHennebelle fc Froman el 2008; H ennebelle fc Tevssier 



1928): 



2008). RAMSES is an extensive code, which can be run in 



many different configurations. We adapt the code in a way 
that it resembles evora (the code we developed and used 
in Paper I) as close as possible. Doing so, we are able to 
compare our findings from one- and three-dimensional sim- 
ulations. In the following we discuss our particular choices 
and discuss differences to the algorithms of evora: 

(i) Similar to evora we use the MUSCL scheme together 
with the MINMOD slope-limiter and the HLLC Riemann- 
solver. To overcome the high-Mach-number problem, RAMSES 
does not use the mechanism described in Sect. B.2 of Paper 
I. Instead it uses a hybrid conservative/primitive scheme 
called pressure fix. In highly divergent fluxes it computes 
the internal energy from the primitive variables. 

(ii) The main advantage of RAMSES (in the context of this 
work) is its ability to adaptively refine the resolution in cer- 
tain regions. We use two complementary refinement strate- 
gies. The first criteria is the mass inside a cell. A cell is 
refined, if the mass inside this cell exceeds eight times the 
mass contained in one cell on the coarsest grid at background 
density. This criteria ensures a higher resolution in denser 
regions. However, it is not able to ensure a adequate reso- 
lution in regions of intermediate density and lower temper- 
ature. Therefore, as a second criteria, a cell is refined if the 
local Jeans length becomes smaller than ten times the size 
of the cell. 

(iii) Due to it's adaptive nature, RAMSES uses a Multi- 
grid Poisson solver to compute the gra vitational potential 
l|Press et af1ll992l ; iKravtsov et aJ1ll997l ). The stopping cri- 
teria parameter of the Poisson solver is set to e = 10~ 6 . 

(iv) The public version of RAMSES offers some very sophis- 
ticated models of c ooling and heat ing, star-formation, and 
feedback processes (|Dubois fc Tevs sier 2008), which are too 
specific and have to many open parameters to be used in 
the context of this study. In order to ensure comparability 
to the evora simulations and to keep the full control over 
the action of the considered physical processes, the cooling- 
module of RAMSES was replaced by the same cooling model 
we use in evora. 

(v) We use the same heuristic constraints on the different 
time steps as in Paper I (see Appendix B therein): C c a = 0.5 
for the CFL-condition and C a — 0.01 for the constraint lim- 
iting the cosmological expansion in one time step. The chem- 
ical subcycling uses also C CO oi =0.1 for the computation of 
the local chemical evolution. 

(vi) In addition to the described changes, a number of 
smaller adjustments have been made in order to use the 
code with the particular physics and initial conditions of 
this study. 



2.3 Jeans length 

In their seminal paper. iTruelove et al.l II 19971 ) demonstrated, 
that in grid based simulation codes which include hydrody- 
namics and self-gravitation perturbations of purely numer- 
ical origin can cause artificial fragmentation. This process 
can be avoided, if the spatial resolution of the code_is suffi- 
cient to resolve the local Jeans length of the gas (|Jeansll 19021. 



Aj = 



Gpto 



n yp 



Gf Ptotal Pga 



(1) 



where c s denotes the speed of sound of the gas. The fraction 
of the local Jeans length and the size of a spatial resolution 
element, called Jeans number J = Aj/Ax, has to be larger 
than one. ITruelove et al.l (|l997l ) sug gested a Jeans numbe r 
of J = 4. In a cosmological setting, ICeverino et al] (|2010h 
found that a Jeans number J — 7 is needed. For our par- 
ticular setup, one-dimensional test simulations show that a 
Jeans number of J = 10 is appropriate. 

Using an AMR code, we can use the Jeans length as 
refinement criterion, increasing the resolution in cold and 
dense regions if necessary. This can be done, however, only 
to a certain level of refinement. In order to be able to pre- 
vent artificial fragmentation on the finest resolution level, 
the standard way is to i ntroduce an artifici a l pressure floor 
llMachacek et al.l l200ll : lAgertz et all 120091 : ICeverino et al] 
120100 . Solving Eq. ([TJ for the pressure, leads to an equa- 
tion for a pressure floor: 



A 2 

Pfloor = Al ptotalPgas 

7T7 



(2) 



If, during the simulation, the pressure in a cell becomes 
lower, the code forces the pressure to this value. Of course, 
this has consequences on the outcome of the simulation: 
Small scale features in the density distribution will be lev- 
eled out. However, in this study we are not interested in 
these small structures. 



3 THREE-DIMENSIONAL COLLAPSE 

3.1 Gravitohydrodynamics 

Similar to Paper I, we begin our considerations with simula- 
tions without the inclusion of radiative cooling and heating 
by the UV background. We start to run the simulation at 
redshift of z = 99. The perturbation on the density field 
is now decomposed into three sinusoidal perturbations each 
along the coordinate axes: 



h A - 



cos (k x x) + A y cos (k y y) + A z cos (k z z) + p) 

(3) 

The corresponding velocity field is obtained using linear per- 
turbation theory: 



fa 



A x sin(k x x) A y sm(k y y) A z sin(fc z z) 



(4) 

We choose the perturbation scale L to be the same for each 
wave, and, again as in our one-dimensional consideration, 
identical to the size of the computational domain in each di- 
rection realizing periodic boundary conditions, at once. The 
choice of that highly symmetrical situation is on purpose 
by reasons initially discussed. In order to compare with our 
one-dimensional simulations we choose identical amplitudes 
along each direction of A — A x = A y = A z , which corre- 
sponds to the choice of identical eigenvalues for the initial 
deformation tensor. The set of initial conditions in the case 
considered here is unlike the situation as for cosmological 
simulations. The cosmological simulations generate for each 



Filaments and sheets of the WHIM 5 




Figure 1. Rendering of the density field at z = for a three-dimensional simulation using L = 4 Mpc. Overdensities in red, underdensities 
in blue. [See the electronic edition of the Journal for a color version of this figure.] 



perturbation scale oc 1/k a set of perturbations at random 
amplitudes distributed in space realizing a given power spec- 
trum of initial density perturbations P(k). Here we consider 
the evolution of one particular realization for a single density 
perturbation on given scale L, only. Therefore, this pertur- 
bation must be given a well-defined amplitude A. This is 
comparable to the perturbation modes with largest length 
scales in the cosmological simulations since these modes are 
also sparsely sampled. However, the value of A affects solely 
the moment of time of the shock formation. The time of 
shock formation is closely related to the moment of time 
when a caustic appears for a one-dimensional collapse of 
collisionless particles provided the length scale and pertur- 
bation amplitude are the same as for the here considered 
collapse of gas. The relation between the perturbation ampli- 
tude and the moment of caustic formation can be obtained 
from the Zeldovich approximation. At first, we set the am- 
plitude equal to A — 0.02 in order to obtain a configuration 
which produces a shock right after the moment 2 = 1 (com- 



patible with the WHIM). In order to investigate the possible 
influence of the value of A, we also perform several simula- 
tions using A = 0.03. These configurations produce shocks 
earlier in cosmic time, and show a more evolved structure 
at Z — 0. Roughly speaking, the output of these two sets 
of simulations is nearly identical but delayed by a redshift 
interval of Az = 1. A more detailed analysis will be given 
below. 

The simulations are performed using a coarse grid of 
256 3 cells and up to 4 levels of refinement. A cell is refined, 
if the mass inside this cell exceeds eight times the mass con- 
tained in one cell on the coarsest grid at background den- 
sity. No Jeans based refinement criteria is applied, yet. We 
include, however, a pressure floor ensuring the Jeans criteria 
(for J = 10) on the finest grid. 

As in the one-dimensional case, the collapse starts with 
an adiabatic contraction of each density perturbation along 
its initial direction producing a two-dimensional sheet. Along 
the line of intersection of the sheets the densities get super- 
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Figure 2. Outcome at redshift z = for a simulation without radiative cooling and heating by the UV background (Panels A-D) and 
with radiative cooling and heating by the UV background (Panels E-H). The simulation uses an initial perturbation length of L = 4 
Mpc. Top Panels: A slice through the filament of 1 Mpc/h X 1 Mpc/h the icy-plane at the position most distant to the halo; Bottom 
panels: A slice through the halo of 1 Mpc/h X 2.5 Mpc/h the xz-plane; Panels A,B,E,F: Density; Panels C,D,G,H: Temperature. 
The color coding of the temperature plot is chosen in way to distinguish photoionized (blue), warm-hot (purple and red) and hot gas 
(yellow). The lines in panel E correspond to the cuts of Fig. [3] The rectangles in panel E and F denote the cuboid volume discussed in 
Sect. [4T4l [See the electronic edition of the Journal for a color version of this figure.] 



imposed, thus producing an elongated filament. The three 
filaments intersect in one point, leading to a further increase 
of density resulting in a halo. The formation of all these 
structures is accompanied by shocks if the initial perturba- 
tion scales are large enough. Since inside the filaments the 
density increases faster than in the sheet (also with respect 
to the one-dimensional collapse), shocks are able to develop 
much faster there. For the same reason, in the sheet far away 
from the filament, the collapse is delayed. The density in the 
halo increases even faster than in the filament, thus the first 
shock emerges here. Outside these structures (halo, filament, 
sheet), the space becomes successively devoid of gas and the 



gas density is lower than the cosmic mean. Thus we will call 
that region widQ 

For illustration, we show a three-dimensional rendering 
of the density field for a perturbation length of L — 4 Mpc 
in Fig.[T] All the described features, the halo, the filaments, 
and the sheets can be noticed. In Fig. [2] A - D, we show the 
density and temperature distribution of the gas in sectional 



2 Throughout the further discussion, we will use the nomencla- 
ture (void, sheet, filament, halo) for the above described struc- 
tures arising in our simulations. We want to emphasize that due 
to the idealized nature of our simulations, these names should not 
be directly attributed to the objects, which are usually labeled by 
these terms in the literature. 
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Figure 3. Profiles through the filament for an initial perturbation amplitude of A = 0.02 and for a series of perturbation scales of 
L = (3, 4, 4.5, 5) Mpc/h and for an amplitude of A = 0.03 for a perturbation scale of L = 4 Mpc/h. Shown are profiles at z = including 
radiative cooling and heating by the uv background, following a line located entirely within the central sheet region (solid lines), and 
following a line coming from the void region and crossing the filament (dashed lines). The direction of these lines is shown in panel E of 
Fig.[2]using dashed lines. For comparison the corresponding one-dimensional pancake profiles from Paper I are also shown (gray lines). 
The profiles are displayed using logarithmic coordinate axes. Upper row: Density profiles; Lower row: Temperature profiles. 



drawings through the simulation box. (Note: the presented 
sectional drawings are projections with depth of one simu- 
lation cell, i.e. we show slices of one cell thickness). The top 
panels show sections perpendicular to the filament's axis (in 
xy-direction), at the position most distant from the halo, 
while the lower panels show sectional drawings along the 
filament (in a;,z-direction) through the halo. The gas distri- 
bution can be distinguished by different phases associated 
with the above described structures: 

• Halo: hot and dense gas confined by strong shocks sep- 
arating the halo gas from the void region. In addition, some 
gas at lower temperature is inflowing along the attached fil- 
aments. 

• Filament: warm-hot gas at I0 4 — I0 5 K with over- 
densities of about 10-100 

• Sheet: cold gas < I0 4 K and low over-densities 8 < 10 

• Void: under-dense gas at even lower temperatures 

Although only gravitohydrodynamics with a polytropic 
cquation-of-state have been assumed so far, streaming fea- 
tures form similar to those which can be noticed in much 
more elaborated simulations of galaxy formation. In our 
case, these cold streams are however very narrow and, near 
the center of the halo, unphysically extended. This happens 
by the finite resolution and the introduced pressure floor. 



3.2 Cooling and heating 

Following the same schedule as in Paper I, we expand our 
simulations by the inclusion of radiative cooling and heat- 
ing by the UV background. Radiative cooling is able to pro- 



duce cool and dense regions. Therefore, in addition to the 
mass based refinement we also apply the discussed refine- 
ment based on the Jeans-length, now. Otherwise we use the 
same parameters as in the previous section. 

In Fig.[2]E - H, we present section slices with respect to 
the outcome of such a simulation for L = 4 Mpc/h similar to 
the non-radiative case discussed above. Several differences 
are evident. The heating by the UV background leads to 
more extended sheets. The filament, forming at their inter- 
section, shows now an extended isothermal core in the cen- 
ter. The filament core has a temperature of ~ 3 x 10 4 K and 
it exhibits very similar properties compared with the core 
region forming in the result of the one-dimensional collapse 
(cp. Paper I). The gas belonging to the cold core streams co- 
herently toward the halo center. During the motion toward 
the halo, the cross section of the stream, i.e. of the core, gets 
narrower and narrower. This cold stream is more extended 
than in the non-radiative case, and penetrates the halo more 
effectively. 



3.3 Pressure floor 

As discussed in Sect. 12.31 a pressure floor has been intro- 
duced to fulfill the Truelove criterion if necessary. This af- 
fects only a small region in the very center of the halo and 
leads to an unphysical high temperature, there. This con- 
cerns also the shape and behavior of the cold streams in the 
vicinity of the halo center. With increasing resolution, the 
necessary pressure floor is lower. As a consequence in higher 
resolved simulations, the streams get even narrower toward 
the halo center. In any case, at these densities our assump- 
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Figure 4. Total halo mass as a function of the initial perturbation 
length for an initial amplitude of A = 0.02 at z = (solid line 
and circles), for A = 0.03 at z = 1 (squares), and for yl = 0.03 
at z = (triangles) . For the computation of the mass, a sphere 
around the halo having a mean density of p/p = 100 is considered. 



tion of an optical thin gas will break down eventually, and 
radiation transport effects would become important. Those 
could work in both directions, expanding the stream by an 
additional pressure support since energy can not longer es- 
cape from the cold stream gas, or they cold constrict it even 
more since UV background radiation gets shielded. 

A turbulent motion in the shocked region which is vis- 
ible in all halo plots is a numerical issue. Though in princi- 
! ile, curved shock s can lead to the occurrence of turbulence 
Paul et al.l l201lh . our highly symmetric initial conditions 
are not expected to lead to a final configuration which ex- 
hibits this kind of motions. Despite the obtained high refine- 
ment level the finite numerical resolution on a rectangular 
grid together with high density gradients in the vicinity of 
the halo center seem to excite these non-physical motions. 
By varying the resolution and the pressure floor, we have 
proven that neither a significant influence onto the dynam- 
ics nor onto the thermal state does exist within the region 
beyond the innermost cells. Nevertheless, these numerical 
shortcomings have to be taken into account if considering 
any results concerning the very central part of the halo. The 
filament regions distant from the halo are not affected, at 
all. 



3.4 Mass of the halo 

The amount of matter contained in the void, the sheets, 
the filaments, and the halo is determined by the initial per- 
turbation scale L. In Fig. [4] we show the mass of the halo 
forming in our simulations. Displayed is the total mass of 
the halo including gas and DM contributions (following our 
assumptions described above). For an initial amplitude of 
A — 0.02 and at a redshift of z — we obtain masses of 
2 x 1O 11 M < M Ha io < 10 13 M Q . Roughly the same behav- 
ior can be found for an higher initial amplitude of A = 0.03 
at an earlier stage of evolution at z = 1. At z — the higher 
amplitude obviously leads to higher halo masses. Since the 
infall of matter onto the halo is still in the linear regime of 



cosmological structure formation, the redshift dependence 
of the halo mass can be well described by the cosmologi- 
cal growth-factor. Since the dependence on redshift and ini- 
tial amplitude obey almost strong proportionality the cor- 
responding curves appear as to be shifted along the mass 
axis. 



4 FILAMENT PHYSICS 

4.1 Filament profiles 

In this section we will discuss the properties of the filament 
distant from the halo. In Fig. [3] we show profiles of density 
and temperature along lines perpendicular to the filament 
and intersecting the filament at the point most distant to 
the halo position. The profiles given by the solid lines are 
obtained by following a line in a;-direction embedded en- 
tirely within the central sheet region, while the profiles given 
by dashed lines are obtained along the diagonal line in in- 
direction, coming from the void region. These lines are lo- 
cated in the projected slice displayed in the upper panels of 
Fig. [2] E, and are indicated by the horizontal line for the 
x-direction, and the diagonal line for the xy-direction, re- 
spectively. Since radiative cooling and heating by the UV 
background introduces an intrinsic physical scale, we dis- 
play a set of simulations with different initial perturbation 
lengths L — (3, 4, 4.5, 5) Mpc/h. The initial perturbation 
amplitude controls the onset of non-linear evolution and 
subsequent shock-formation in the simulations. Therefore, 
we have also included a simulation with a higher initial am- 
plitude of A — 0.03 into our considerations. Thus we obtain 
a more evolved structure at z = 0. 

The profiles show similar properties as the one- 
dimensional pancakes of Paper I. They are, however, more 
complex and the reached densities are by an order of magni- 
tude higher. Provided an initial perturbation length of ~ 4 
Mpc/h a shock is able to form. Then, the obtained profile 
can be separated into an isothermal core of 10-100 kpc/h 
and a surrounding shock-heated region. The spatial exten- 
sion of the core is again determined by the perturbation 
scale L. Larger initial scales L lead to smaller cores. The 
relation, however, appears to be much steeper compared to 
the one-dimensional case. This results in a smaller possible 
scale range for L where a core can exist. 

The profiles across the filament obtained for the dif- 
ferent directions (one along the sheet and the other start- 
ing from the void region), show no differences for the inner 
(isothermal) part. For the boundary region of the core how- 
ever, the density profile along the sheet shows an almost 
smooth transition from the sheet to the filament region, 
while for the profile following the direction from the void, a 
step-like increase can be noticed. In the first case, the cor- 
responding temperature profile exhibits a rather continuous 
increase toward the filament region, but for the second case, 
a clear shock is noticeable. 



4.2 Scaling relations 

In Fig. [S] we show the dependencies on the initial pertur- 
bation scale L for the temperature T c and for the number 
density of hydrogen nuc in the center of the filament. Shown 



Filaments and sheets of the WHIM 9 



5 


;'l 


\ ' 


A 




-3 


-'i 


, , ,,_, , ^ 


2 4 

o 


-A 
■ 


■ \ 




CO 

E 
.u. 

c 

Ol 

o 


-4 
-5 
-6 


A 
.■ 


■ 

A ■ 

A r 

.r ; 


3 










-7 






1 


10 




1 


10 




log(L) [Mpc/h] 








log(L) [Mpc/h] 


-14 










5 








A / 


C 


" ' 


r»- 








E 
o 
c -16 

-Q-17 

Q. 

15) 
o -18 


~A 


A / 

A ■/ 
A M 
A M 




o 
H 

Ol 

o 


4 


: ■ 


^^\ A^ A 

■ A - 


-19 










3 






1 


10 




7 


-6 -5 -4 






log(L) [Mpc/h] 










log(n Hc ) [cm" 3 ] 



Figure 5. Central temperature T c and the central hydrogen 
number density njj c of the filament for an initial amplitude of 
A = 0.02 shown at z = {black squares) and for an initial am- 
plitude of A = 0.03 shown at z = 1 (open triangles). Top Pan- 
els: Dependence on the perturbation length scale L [Panel A: 
Central temperature; Panel B: Central hydrogen number den- 
sity; Panel C: Central thermal pressure). The corresponding one- 
dimensional results from Fig. 6 of Paper I along with the derived 
one-dimensional scaling relations are given for comparison {grey 
points/line). The three-dimensional scaling relations for the fila- 
ment, derived in this paper, are indicated by black lines. Panel D: 
Phase space diagram, i.e., the dependence of the central tempera- 
ture versus the central number density {black points) for different 
perturbation length scales L. The solid black line denotes the 
equilibrium temperature T Gq at redshift z = and the dashed 
line denotes the equilibrium temperature T oq at redshift 2 = 1. 
Again, the corresponding one-dimensional results are given for 
comparison {grey points). 



are results for two sets of simulations: our reference model 
with an initial amplitude of A — 0.02, and, for comparison, 
a simulation with a higher initial amplitude of A — 0.03. 
The latter shows faster evolution with time. We therefore 
compare the output of the simulation employing A — 0.02 
at z = to the A = 0.03 simulations at z = 1. We compare 
the obtained relations for the three-dimensional case with 
those obtained for the one-dimensional case as discussed in 
Sect. 5 of Paper I. Filaments on scales L < 4 Mpc/h do 
not exhibit any shock and their general behavior is simi- 
lar to what we found for the one-dimensional pancakes on 
scales less than 2 Mpc/h. For perturbations evolving over 
larger scales, shocks appear and significant differences be- 
tween the three-dimensional and one-dimensional cases can 
be noticed. In both cases, an increasing scale length L leads 
to larger central densities, but with a much stronger depen- 
dence for the three-dimensional case. The stronger increase 
in density is accompanied by a stronger decline in tempera- 
ture. This can be noticed in panel A and B of Fig. \5\ In the 
three-dimensional case the T C (L) and nn c (L) relations are to 
some extent compressed with respect to the one-dimensional 



case. Roughly speaking, the one-dimensional scaling relation 
in the interval 2 Mpc/h < L < 24 Mpc/h corresponds in the 
three-dimensional case to a much narrower interval 4 Mpc/h 
< L < 5.5 Mpc/h. A corresponding behavior can be found in 
the phase-space diagram given in panel D of Fig. [5] There, 
we plot the central temperature versus the central hydrogen 
number density. The behavior of the T'c(w-Hc) relations differs 
only marginally for the one- and the three-dimensional cases. 
Therefore, for the given interval of L, we can apply the same 
arguments for the explanation of the observed scale depen- 
dencies as given in Paper I. In particular, the temperature 
of the core is determined by the balance of radiative cooling 
and heating by the UV background. We call this equilibrium 
temperature T eq . For the simulations employing A — 0.03 we 
observe the same behavior at z — 1. At this point in cosmic 
time, however, our model for the UV background produces 
a flux which is by an order of magnitude stronger, result- 
ing in a different equilibrium temperature. Nevertheless, the 
dependence on the perturbation length remains unchanged. 
Since the physical densities are plotted the higher values for 
density at z — 1 (A — 0.03) is due to accounting for the dif- 
ferent scale factor at higher redshift. For small L, as in the 
one-dimensional case, the reached densities are to low for 
effective cooling at any time. Thus, after the initial heating 
by the UV background, the thermal state in the center is 
determined by the interplay between adiabatic contraction 
of the gas and the cosmic expansion, only. 

For an one-dimensional collapse, it can be shown 
that the temperature of the shocked gas scales always as 
T oc L 2 (cp. with the re sults of our paper I; see also 
IjSunvaev fc Zeldovichlll972l ). In the result, the thermal pres- 
sure of the infalling gas, surrounding the filament, scales 
also oc L 2 . This is also true for the thermal pressure inside 
the sheets, since these are the result of a one-dimensional 
collapse process. The thermal pressure in the core of the 
filament, however, is in equilibrium with the the dynamic 
pressure 



Pdyn, sheet 0*. "He, sheet ^c, sheet 



(5) 



of the gas compressed into the sheet and streaming onto 
the filament's core. For the velocity component toward the 
filament u c , sheet, one obtains a scaling behavior of it Cj shcct oc 
L from basic hydrodynamical considerations. As shown in 
Paper I, the scaling behavior of riHc,sheet can be derived from 
the density dependence of T eq . For densities of p/p — 5 — 
15, appropriate for what we find in the sheets, we obtain 
T cq oc rig°' . As already discussed, the gas pressure outside 
the sheet scales according to p c , sheet oc T cq nHc, sheet oc L 2 . 
Following the same line of arguments as in paper I, this leads 
for the sheets to riHc, sheet oc L 25 . The state of the core of 
the filament is determined by the dynamical pressure of the 
sheet, acting on the core of the filament. It scales therefore 
like 



r 4 
Pdyn, sheet OC Lj 



(6) 



since u c , sheet oc L. This relation determines the pressure 
inside the core of the filament, which therefore shows a much 
stronger scale dependence as for the core layer of the one- 
dimensional pancake. This can be noticed in panel C of Fig. 
[S] Since the temperature inside the core is determined by the 
equilibrium temperature, the density is the only variable to 
adjust a nearly hydrostatic equilibrium. Hence, the relation 
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between riuc(L) for the core of the filament is much steeper 
than for the one-dimensional pancake. In order to obtain 
scaling relations for the filament, we can apply the same 
consideration as before. Since the cores of the filaments are 
denser than the sheets (p/p = 100 — 500), we have to use a 
slightly different slope for the equilibrium temperature T cq oc 
n Hc 21 - Using Euler's equation for the estimate of the central 
pressure p c oc riHc4>, and Poisson's equation <f> oc whcAq, 
we obtain the relation p c oc X^n^, where A is the spatial 
extend of the core. Using the scaling relation for the pressure 
inside the filament's core p c oc L 4 ' 5 derived above, we obtain 



TiHc 



T c oc U 
A oc L 



(7) 
(8) 
(9) 



These analytically derived dependencies approve the out- 
come of our simulations shown in Fig. [5] This justifies our as- 
sumption about a thermo-dynamical equilibrium state of the 
filament's core. The relations obtained for the 3-dimensional 
filament completely differ from the results obtained for the 
one-dimensional pancake. Though the n — T relation still 
holds for the equilibrium state of the core where radiative 
cooling and heating due to the UV-background flux balance 
each other we get partly much stronger dependencies on 
the initial perturbations length scale. This is a direct conse- 
quence of the 3-dimensional gas flow structure. In addition, 
the analytical findings allow for extrapolating to values at 
larger L where the necessary numerical resolution cannot be 
reached anymore. 

As direct consequence of the much stronger dependence 
Xo{L), a core of reasonable spatial size can only exist for a 
significant smaller range of perturbation scales. If we use the 
L — 4 Mpc/h case with A = 60 kpc/h as a reference, we 
obtain a scale of 



L= 14 Mpc/h 



(10) 



where the core of the filament becomes smaller than 1 kpc/h. 
Our above made assumption about the tight coupling of 
baryons and DM might lead to local deviations from the de- 
rived scaling relations within the core region. As discussed 
in Sec. [2] admitting the full initial perturbation spectrum, 
clustering of the DM on small scales may happen unlimited. 
These DM clumps may trigger local instabilities of the gas 
to collapse. These collapsed objects are expected to reach hy- 
drodynamical equilibrium with their ambient medium after 
some time. The thermodynamical conditions of the ambi- 
ent medium (core region of the filament) however are deter- 
mined by the above described processes. Thus, it can be ex- 
pected that the given scaling relations are valid throughout 
the filament's core region except for the very local collapsed 
regions. With other words, averaging over scales larger than 
the possible gravitational instabilities the derived scaling re- 
lations should hold. To follow the evolution of the mentioned 
instabilities in detail is still beyond the numerical feasibili- 
ties. 



4.3 Formation process 

In order to discuss the detailed formation process of the fil- 
ament we use the simulation employing a higher amplitude 
of A — 0.03. For this initial amplitude, the filament reaches 



the stage of non-linear evolution and shock-formation before 
2 = 1. Thus we are able to follow the formation of the fila- 
ment for a longer time and up to a later stage of evolution. 
In Fig. [6l we show the temperature distribution in slices 
analogously to the upper right panel of Fig. [2] at different 
moments in time. The left-hand panel shows the evolution- 
ary stage at the redshift z = 1 right after the shock has 
formed. The subsequent evolution is shown on the middle 
panel at z = 0.667, and the right-hand panel at z — 0. 

Before shock formation, matter accretes onto the fil- 
ament almost equal concentrically but slightly enhanced 
along the planes given by the sheets where the density is 
higher. Density and temperature in the filament increase 
adiabatically, until a shock begins to form. Because of the 
much lower sound speed in the void region, a shock is form- 
ing first at the interface between the forming filament and 
the void. Further accretion can be distinguished by two 
phases: 

• Material accreting from the void, which is shock-heated 
to warm- hot temperatures (T > 10 5 ). 

• Material accreting along the sheet, which is heated by 
the transformation of its gravitational energy if the gas is 
moving toward the center, only. 

This corresponds to our findings in Sect. 14. l\ where the pro- 
file exhibits a shock at the interface to the void region but 
not in direction to the attached sheet. Noticeable, the ve- 
locity field is a good tracer for the shock, declining rapidly 
behind the discontinuity. 

The isothermal-core of the filament is not influenced by 
this further accretion. The material accreted from the sheet 
is diverted toward the void resulting in an effective outflow 
from the filament. This leads to a cold diffuse clover-leaf- 
shaped region around the filament. A similar behavior is 
found in the one-dimensional pancakes, where the subse- 
quent accretion also does not influence the central distri- 
bution, but rather accumulates gas around the core. The 
observed outflow however is only possible in two or more 
dimensions. Summarizing, the gas inside the filament can 
be distinguished into three-phases: a warm-hot shocked gas 
layer, an isothermal core, and a cold diffuse region in be- 
tween. After some moment of time, the core gets "shielded" 
against further infall of gas and remains thus a limited reser- 
voir of cold gas streaming towards the halo. Besides the de- 
scribed reflection of infalling gas, the hot gas in the immedi- 
ate vicinity of the cold core is not able to cool anymore due 
to the large cooling time. 

4.4 Time evolution of the core 

Next, we are considering the time evolution of the amount 
of gas which resides in the core of the filament. To that 
purpose, we extract a cubical volume of size 0.6 Mpc/h x 
0.6 Mpc/h x 1.0 Mpc/h embedding a representative part 
of the filament region. This cuboid is centered in x- and in- 
direction. In z-direction it covers the range from the edge of 
the box up to half of the distance to the halo center. The 
rectangular projections of the cuboid are shown in Fig. [2] 
E and F (dashed lines). For the chosen volume we compute 
the fractions of the gas at different states. W e first construct 
a ph ase-space diagram (T versus p/p, cp. Valage as et al.l 
120021 ). weighted by the mass of the individual cells. Then, 
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Figure 6. Time evolution of the temperature distribution of the forming filament. Shown are snapshots at expansion factors of a = 
(0.6,0.8,1-0) corresponding to redshifts of z = (1.0,0.6667,0). The range of the plot is similar to Fig. [2]E. It uses, however, a linear scale. 
The vectors visualize the gas velocity in the plane of the slice. The white line denotes an over-density of 10 1 ' 5 . [See the electronic edition 
of the Journal for a color version of this figure.] 



we separate the gas contributions according to the differ- 
ent phase states. We use a density of p/p = 10 1 ' 5 as the 
threshold marker in order to distinguish between the true 
core gas and the cold and diffuse outflow region. In Fig. [7J 
we display the time evolution for the dense gas forming the 
isothermal core as well as the less dense envelope around it 
(computed from subsequent simulation snapshots). The evo- 
lution curves are displayed for simulations with L = (4, 5) 
Mpc/h for an initial amplitude of A — 0.02, as well as for the 
simulation with L — 4 Mpc/h at the higher initial amplitude 
A = 0.03. 

For the simulations employing the lower initial ampli- 
tude, the fraction of dense gas, displayed in the upper panel, 
rises monotonically but with decreasing slope till a = 1. For 
larger L more cold gas accumulates in the core, which is 
in agreement to our findings in Sect. 14.21 The higher initial 
amplitude of A = 0.03 allows to follow up the evolution up 
to later stages. From Fig. [Jjone can read, that the use of a 
higher initial amplitude is shifting the plot for L = 4 Mpc/h 
to the left along the time axis. The plot shows that the frac- 
tion of the cold and dense gas decreases at later moments of 
time. 

For the fraction of gas in the gaseous envelope, displayed 
in the lower panel of Fig. [7] we observe an inversed behavior 
in respect to the mass fraction of the gas inside the core. The 
part of the envelope, which has a temperature of T > 10 4 ' 5 
and can therefore be considered warm- hot, however, follows 
roughly the evolution of the core. Deviations are caused by 
the complicated temperature distribution around the core. 
When the cold core emerges in the center of the forming fila- 
ment, the fraction of the warm-hot enveloping gas is decreas- 
ing. Later, when the core is formed and is shielded against 
further infall, the mass fraction of the ambient warm-hot gas 
is increasing, again. At later stages, a decrease of the warm- 
hot fraction can be noticed. This is caused by the adiabatic 
cooling due to the cosmological expansion. 

Hence, from a given point in time, the filament evolu- 
tion leads to a continuous decrease of the mass fraction of 



the isothermal core. This stage of evolution, however, is only 
reached untill z — for filaments forming from perturba- 
tions with a high enough initial amplitude. The process can 
be explained by the process described in Sect. 14.31 Starting 
from the denser regions next to the halo, the shock encloses 
more and more of the filament, until no cold gas can accrete 
on it from the void. The remaining reservoirs of cold gas are 
the sheets. However, this gas does not reach the core any 
more, but it is redistributed around it. Since, with time, the 
cold gas in the core is drained toward the halo, the amount 
of cold dense gas in the filament is decreasing. 



5 ACCRETION ONTO THE HALO 

5.1 Formation process 

In this section, we consider the formation and the evolution 
of the collapsed structure forming in the center of our simu- 
lation box which can be assigned to a gaseous halo. This halo 
forms at the intersection of the three filaments. Concerning 
the halo, our assumption about the identical distributions 
of gas and DM must be considered with caution. During its 
evolution, the gas halo could undergo a significant concen- 
tration in comparison to the DM distribution. In particular, 
this may happen due to the cooling processes. Shock heat- 
ing, the introduction of an artificial pressure floor, and heat- 
ing by the UV background acts against that concentration, 
however. On average, the halo gas enclosed by the shocks 
follows a nearly isothermal distribution. Thus, the two dis- 
tributions (DM and gas) should not differ significantly from 
each other. Concordantly, the potential in our simulations 
corrected for the gravitational action of the DM is still a 
sufficient good approximation and there is no necessity to 
fully include the DM into our simulations. 

With these preconditions in mind, the halo region of 
our simulation box offers nevertheless an excellent testing 
ground to study the hydrodynamical effects influencing the 
accretion over time (exclusive of more complicated baryon 
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Figure 8. Time evolution of the temperature distribution of the halo. Shown are snapshots at expansion factors of a = (0.4,0.5,0.6) 
corresponding to redshifts of z = (1.5,1.0,0.6667). Again, the figure uses a linear scale. The vectors visualize the gas velocity in the plane 
of the slice. [See the electronic edition of the Journal for a color version of this figure.] 



physics, like star formation, fragmentation, turbulent mo- 
tions etc.). In particular, the supply of the inner halo re- 
gions with cold gas can be studied in this way. The forma- 
tion processes for the halo exhibits some differences to the 
formation process for a filament being described in Sect. 14.31 
The density of the forming halo is much higher than inside 
the filament, and as a result, much higher temperatures are 
reached. This is caused the by strong heating due to adi- 
abatic compression and by stronger shocks, caused by the 
deeper potential. 

Fig. [8] shows the time evolution of the temperature dis- 
tribution in the halo region. Initially, a hot and dense clump 
forms by adiabatic compression at the intersection of the 
filaments. When the potential well of this clump is deep 
enough, a shock forms at the interface toward the void re- 
gion. This shock is moving outwards with time and affects 
also the adjacent filament regions. The high pressure of the 
surrounding halo medium leads to a constriction of the in- 
flowing filament gas. Therefore, the accretion onto the halo 
can be distinguished by two modes: 

• The gas continuously infalling from the void region onto 
the halo is shock-heated to high temperatures T > 10 6 K. 

• The colder gas, forming the isothermal core of the fila- 
ment, can propagate deep into the halo supplying the halo 
with cold gas. 

The shock around the halo forms earlier than the shock con- 
fining the filaments. Only at later times, the subsequent gas 
accreting onto the latter produces shocks. Those shocks con- 
fine the cold isothermal cores of the filaments. 

As already mentioned, the described mechanism shows 
strong similarities to the cold stream scenario. We can now 
identify the cold streams with the isothermal cores in our 
simulated pancakes and filaments. The formation and evo- 
lution history obtained by our simulations leads immediately 
to an important conclusion: The cold streams have not to 
penetrate the shock around the halo, but are forming and 
existing at the very beginning of the evolution. They are 
attached to the halo even before any shock fronts around 



the halo have formed. The later forming shocks are influ- 
encing only the shape of the already existing streams, i.e. 
the cross section of the attached filaments and of their cold 
cores. This compression of the filament gas streaming into 
the halo was already shown in fully cosmological simulations, 
e.g. lDekel et al.1 (|2009T i . iKeres et ail j|2009bT ). 



5.2 Quenching by thermal conduction 

According to the above findings, the cross sections of the 
filament cores consisting of cold and dense gas are decreas- 
ing with increasing initial perturbation scale L. The core 
size undergoes contraction along the filament's axis toward 
the gaseous halo and its gas content is time-dependent. In 
Figure [9] we show density and temperature profiles across 
the filament core (in ^-direction, similar to the solid lines in 
Fig. [3} for L — 4 Mpc/h. The profiles are given at differ- 
ent locations near to the halo, where the shocked halo gas 
starts to influence the shape of the cold streams along the 
filament cores. It can be clearly noticed that with decreasing 
distance to the halo center the temperature gradient at the 
boundary of the isothermal core is significantly increasing 
and the cross section of the core is decreasing. High temper- 
ature gradients, however, might be accompanied by efficient 
thermal conduction. As discussed in Paper I, under certain 
conditions this may lead to a partial or entire evaporation 
of an initially cold region. An order-of-magnitude estimate 
gives for the case that thermal conduction may overcome 
the cooling, i.e. —V ■ j > A, where j denotes the thermal flux 
and A the cooling function: 



n e A T < 0.74 x 10 



MM 
Vio 4 kJ 



-*- shock 

10 6 K 



(11) 



where n e is the electron number density, At is the character- 
istic scale for the temperature gradient, T c is a temperature 
of the core of the stream, and Tshock is the temperature of the 
halo confining the stream. Under the prevailing thermal con- 
ditions and due to the UV background the gas is almost fully 
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Figure 7. Time evolution of the mass fraction of the gas which re- 
sides in the isothermal core and has an over-density of p/p > 10 1,5 
(upper panel), and for the gas outside the core having an over- 
density of p/p < 10 (lower panel, upper curves) for the cuboid 
volume around the filament described in the text. Also displayed 
is the evolution of the mass fraction of the gas at temperatures 
T > 10 4 ' 5 K and p/p < 10 1,5 (lower panel, lower curves), which 
characterize the evolution of the hot gas envelope around the 
core. Shown arc the curves for the simulations using L = 4 Mpc 
(dashed line), L = 5 Mpc (dotted line), and L = 4 Mpc, but with 
an initial amplitude of A = 0.03, thus forming shock earlier (solid 
line). 



ionized n e w nu- Saturation does not take place. The pres- 
ence of a static magnetic field may considerably lower the ef- 
ficiency of Spitzer's relation. However it is rather likely that 
the magnetic fields are not static but generated by turbu- 
lent chaotic motions. In that case, the heat conductivity can 
approach nearly Spitzer's c oefficient (jNaravan fc Medvedevl 
l200ll ; IVoigt fc Fabianll2004r i. Thus Spitzer's relation for ther- 
mal conductivity is a good approximation. In addition we 
use quantities in order to estimate the lower limit of a pos- 
sible effect. 

The contraction of the cold stream inside the halo is 
caused by the increased pressure, which is a result of the 
collapse process. At the outer boundary of the halo (given 
by the outwards moving shock) the spatial extend of the cold 
stream is still determined by the scaling relations obtained 
in Sect. 14.21 The temperature of the shock of the halo, how- 
ever, is caused by shock-heating alone. It can be therefore 
estimated from the transformation of the kinetic energy of 
the gas into thermal energy, when infalling onto the halo. 



Figure 9. Temperature profiles through the cold stream in close 
vicinity to the halo for an initial perturbation amplitude of 
A = 0.02 and a perturbation scale of L = 4 Mpc/h. Shown are 
profiles at z = including radiative cooling and heating by the 
UV background. Like in Fig. \3\ (solid lines), we follow a line lo- 
cated entirely within the central sheet region. The profiles are 
displayed using logarithmic coordinate axes. Displayed are pro- 
files with decreasing distance to the center of the halo: 340 kpc/h 
(solid line). 240 kpc/h (dashed line), 200 kpc/h (dottet line), and 
180 kpc/h (dashed- dotted line). 



This leads us to 



fc B T s , 



mn 



GpL 2 



T sl 



10 6 K 



0.62 



(- 



-V 



V4Mpc/hy 



(12) 

where mn denotes the mass of the hydrogen atom, fee the 
Boltzmann constant, and p the background density. Thermal 
conduction becomes effective if the transition length At will 
be of comparable size as the cold stream's cross section A 
or larger, i.e. A < At- From Eq. (|11|) (transformed into 
dimensionless units) we obtain 



( fine \ 
Vlcm-3,< 



10 4 K 



Is, 



( A 



io 6 Ky V 1M P c / h 



> 2.3x10" 
(13) 



In Sect. l4~2l we obtained nuc oc L 5 ' 7 , T c oc L , and A oc 
IT 3 ' 57 . From Fig.[9]we obtain A Q « 0.025 Mpc/h, T c « 2.5 x 
10 4 , and (p/p) c » 200 (corresponding to n Hc = 4.93 x 10~ 5 
cm -3 ) at L = 4 Mpc/h. With these values, as well as Eq. 



(|12[) . above inequality becomes: 
\4Mpc/hJ 



< 0.52 



(14) 
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When solving for L, we obtain a critical scale of 

L crit « 6 Mpc/h (15) 

from which on thermal conduction may lead to an evapora- 
tion of the cold stream. This corresponds to a halo mass of 
M Ha io « 1O 13 M . 

The quantities T s hock, u c , sheet, etc., which we consid- 
ered in Sect. 14.21 depend on the redshift at w hich the shock 
formes (see also ISunvaev fc Zeldovichl ll972T ) . Taking these 
dependencies into account we obtain for A oc (1 + z s ) _1 
and for At oc (1 + z s ) -5 ' 4 . Here z s denotes the redshift of 
shock formation. If we include these dependencies into our 
estimate for the critical scale Lcrit, we obtain 



L c 



(1 + Zs 



(16) 



This dependence on the shock formation time is rather mod- 
erate and makes the estimate for L cr it relatively robust. 

Inside the halo, where the temperature gradients are 
even higher, thermal conduction becomes effective already 
at L < L cr it, depending on the distance to the center of the 
halo. This should lead to a decrease of the penetration depth 
of the cold streams. Because of the drastic increase of compu- 
tational complexity, we did not include thermal conduction 
into our simulations self-consistently. However, the above 
considerations together with the computational results in 
Paper I show, that under the existing conditions thermal 
conduction leads to an evaporation of the cold filament core. 
Concordantly, at initial perturbation lengths L > L cr it, the 
formed cold streams are not longer able to penetrate deep 
into the halo, but evaporate if the halo forms a shock and 
starts to constrict the attached filament. Though the an- 
alytical considerations lead to a quite robust estimate for 
the critical L cr it the time-dependent evolution can only be 
obtained by future simulations. 



6 SUMMARY AND CONCLUSIONS 

In the present paper, we extend our investigations related to 
the formation of sheets in the intergalactic medium of the 
Universe as described in Paper I to a fully three-dimensional 
gas distribution. Again, we consider a highly symmetrical 
setup of initial and boundary conditions. We follow the evo- 
lution of one single perturbation mode sized according to the 
dimensions of the simulation box. In Paper I we have shown, 
that under conditions of shock heating the influence of al- 
most all small-scale modes will be erased during the evolu- 
tion of the largest mode provided the gas will be shocked on 
that scale. Taking into account only the largest perturbation 
mode, we are confronted with a problem similar to the cos- 
mic variance problem in any cosmological simulation: The 
perturbation amplitudes related to the largest modes are 
sparsely sampled or even represented by a single value, only. 
Thus we chose an amplitude of the considered perturbation, 
which is not far from the cosmic power spectrum. We con- 
sider, however, different values for the initial amplitude, in 
order to show possible differences in the evolution of the per- 
turbations. As in Paper I, we successively include different 
physical processes such as radiative cooling and heating due 
to the UV background into our simulations. In comparison, 
the outcome of our simulations, a three-dimensional con- 
figuration of sheets, filaments, and a gaseous halo, is much 



more complex. In particular, the higher degree of freedom 
leads to a complicated velocity field and well pronounced 
gas flows toward the maximum of gravitational potential. 
The main conclusions we derive from our three-dimensional 
simulations are as follows: 

(i) Formation and structure of filaments: We are 

able to study the density and temperature distribution as 
well as the velocity field of the evolving filaments in great 
detail. Though the resolution of our simulation in the halo 
region is comparable with the obtained resolution in cos- 
mological hydro simulations, the acquired resolution for the 
filaments (at low to intermediate gas densities) is in our 
simulations much higher. For these particular environment, 
our special setup, combined with the additional refinement 
according to the Jeans criterion, enables us to reach a reso- 
lution, which lies about three orders of magnitude in mass 
above the state-of-the-art IGM simulations (see also the 
footnote in the introduction). We get much more insight 
into the detailed structure, in particular the velocity field, 
of filaments of a scale up to 5 Mpc/h. In simulations with a 
lower spatial resolution these details would not be accessible. 
In particular, regions of high density and low temperature 
would be smoothed out by the pressure floor or a compara- 
ble mechanism to prevent a violation of the Truelove crite- 
rion. Less resolution makes necessary a higher pressure floor. 
Therefore those simulations become more and more similar 
to pure adiabatic simulations without coo ling (a similar con - 
clusion was obtained in considerations bv lKang et al.ll2005h . 
For these reasons, the investigation given here might be con- 
sidered as complementary to the existing hydrodynamical 
cosmological simulations. 

The gas streaming onto the filament forms a shock if the 
initial perturbation length exceeds L = 4 Mpc/h. The shock 
does not form everywhere at once, but develops first toward 
underdense regions, and successively covers the whole fila- 
ment. In its center, the shock confined filament exhibits an 
isothermal core at temperatures of T ~ 2 x 10 4 K. As for the 
one-dimensional case, sheets are forming, but they are not 
longer confined by shocks. The streaming motions toward 
the maximum of the gravitational potential cause a delayed 
collapse behavior and suppress shock-formation. 

The filament is provided with new gas essentially in two 
ways. On the one hand, gas accretes directly onto the fila- 
ment (from the void region). This gas experiences a strong 
shock, which temperature is essentially given by the initial 
perturbation length. On the other hand, gas is flowing along 
the sheets towards the filaments. The shocks formed by this 
gas are closer to the filament's center and less pronounced. 
The thermodynamical state of the gas at the center of the 
filament is mainly determined by this gas. 

The time of shock formation is governed by the value of 
the initial perturbation amplitude. A higher amplitude leads 
to an earlier shock. 

(ii) The state and the evolution of the inner core 
of the filament: One of the main characteristics of the 
filament is the state of the gas at its very center. Here, it 
forms a well confined core along the filament's axis consist- 
ing of cold and dense gas. The temperature of this gas is 
determined by the balance between heating due to the UV 
background and radiative cooling. In the early stage of the 
formation of the filament, due to the initial compression of 
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the gas, radiative processes are very effective, i.e. cooling 
and heating times are short. Then, an equilibrium state be- 
tween the (time-dependent) UV flux and cooling processes 
can be reached. During the further evolution, the time evo- 
lution of the UV flux and the overall density dilution due 
to the cosmic expansion moves this equilibrium state along 
the phase trajectory T c (nuc)- The quantities related to that 
core (nHc,Tc,Ac, etc.) exhibit a very strong dependence on 
the initial perturbation length L. For the filaments in the 
three-dimensional case this dependence is much stronger and 
the physical mechanisms are different compared to the one- 
dimensional case. For the gas accreting directly (not via the 
sheets) the temperature scales again oc L 2 . This is a gen- 
eral outcome, and can be derived from the hydrodynamical 
equation s as long as dissipative pro cesses can be neglected 
(see also lSunvaev fc Zeldovichlll972T ). 

It can be noticed from the formation process of the inner 
region of a filament that the state of the core is determined 
by the gas inflowing along the sheets. The dominant gas 
pressure is in this case not longer determined by the gas tem- 
perature in the shock surrounding the filament, but by the 
dynamical pressure of the gas streaming in via the sheets. 
Taking that into account, we are able to derive analytically 
the scaling relations confirming the results obtained from 
our simulations. The initial amplitude does not enter these 
relations. It determines solely the moment of time when the 
shocks and the cold core form. However, since the UV flux jo 
is time dependent and defines the equilibrium state at any 
moment, the moment when the shock forms also defines the 
starting point of the equilibrium track along T c (nuc, jo)- The 
obtained scaling relations allow for a reliable extrapolation 
of the results beyond L = 6 Mpc/h where the obtainable 
resolution becomes unsufficient. In particular, the strong L- 
dependence of the core size A c leads to an extremely con- 
stricted core for large L > 14 Mpc/h at redshifts z ~ 0. 

(iii) The cold streams: The velocity of the gas, which 
belongs to the cold core is aligned along the local filament 
axis and is directed toward the center of the collapsing 
halo. Concordantly, these streams provide the inner halo 
region with cold and dense gas. The overall properties of 
these streams correspond to the cold streams observed in 
hydrodynamical simulat ions of galaxy formation (see e.g. 
iNagai fc Kravtsovl|2003f ). Hence, if one attributes the halo 
to a galaxy, this gas would significantly impact on the star 
formation history of this galaxy. We find these streams to be 
an essential characteristic of our configuration, and they ex- 
ist from the very beginning of the collapse process. Being a 
part of the core of the filament, the stream follows the same 
scaling relations. Therefore, the stream becomes narrower 
and denser with increasing L. The influence of the outward 
moving accretion shock of the halo results in a further con- 
striction of the stream. This effect becomes stronger with 
decreasing distance to the halo center. 

The total amount of gas which can be accreted via cold 
streams, is determined by the same processes which control 
the core of the filaments. The main constituent in this con- 
text, is the interplay between radiative cooling and heating 
by the UV background. The combination of these processes 
with the particular hydrodynamic conditions, given by the 
characteristic scales of the system, result in the scaling rela- 
tions discussed above. Besides the radiative cooling and the 
heating by the UV background, additional processes acting 



on the energetic budget of the gas, e.g. radiative feedback 
from the forming galaxy, may have an impact on the spatial 
extend of the stream as well as on the amount of transported 
gas. 

After the filament has formed, the amount of mass inside 
the filaments core is decreasing with time. The gas infalling 
from the void region onto the filament is shocked and can 
therefore not reach the core region, anymore. The cooling 
time for the gas at high temperature and low density has 
become too large. At some stage of evolution, also the gas 
streaming along the sheets does not reach the filament's core 
region. This happens when the dynamic pressure of the gas 
flow along the sheets becomes equal or less the pressure of 
the cold and dense gas of the filament's core. Therefore, 
the core is shielded from further gas supply. In the result, 
the total reservoir of cold gas available for star formation 
processes form cold stream accretion is limited to a well 
defined fraction. This gas fraction is determined soon after 
the formation of the filaments. With time, the core gas is 
draining off toward the gaseous halo and the mass content 
of the filament cores is expected to decrease continuously. 
This is a consequence of structure formation only. 

Taking into account the full initial perturbation spectrum 
for the DM which was not considered here might alter some 
results. In particular, the gas content in the cold streams 
available for gas supply of the halo may considerably de- 
crease due to the formation of collapsed objects within the 
filament's core region. As in all cosmological simulations, 
the large-scale power is cut due to the fixed size of the sim- 
ulation box. This may have impact on the time evolution. 
Taking into account the action of the large-scale modes of 
the perturbation spectrum beyond the actual box size, tidal 
fields, angular momentum transfer etc. will modify the ve- 
locity field and in particular the geometrical shape of the 
cold streams. Therefore, the gas flow into the halo will also 
change with time. 

(iv) Termination of cold streams: The enhanced pres- 
sure of the shocked gas inside the halo leads to a subsequent 
constriction of the cold streams reaching the inner halo re- 
gion. This effect increases with decreasing distance to the 
hot and dense center of the halo and is caused by the adia- 
batic contraction of the gas. It needs the additional influence 
of thermal conduction to prevent the cold gas from stream- 
ing till almost the center of the gas halo. At the outer vicinity 
of the halo, large temperature gradients occur. These gra- 
dients become even stronger towards the center of the halo. 
As a result, thermal conduction is effective, and the heat- 
ing by the thermal flux toward the cold stream region is 
able to exceed the radiative cooling there. Using the scaling 
relations for the core of the filament, we obtain that for ini- 
tial perturbation lengths L > 6 Mpc/h, which corresponds 
to a halo mass of about M^aio ~ 10 13 Mq, the cold stream 
evaporates at the boundary of the accretion shock and is 
therefore not able to exist inside the halo, anymore. For in- 
termediate perturbation lengths 2 Mpc/h < L < 6 Mpc/h a 
partly evaporation of the streams, caused by constriction of 
the stream inside the halo, may occur. This process should 
start near the central region, first. Thus, a termination of 
cold gas supply may happen for smaller scales as well. The 
corresponding mass range for which the cold streams can 
reach the inner halo is 3 x 10 n M Q < Af ha i < 10 13 M G . 

These results are calculated for z ~ 0. If the initial 
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perturbation amplitude is higher, the shock formation and 
the collapse of the halo start considerably earlier, i.e. at 
higher redshifts. As a consequence, the physical density 
in the filament's core is higher. The characteristic length 
for thermal conduction At decreases with increasing red- 
shift stronger than the core size Ao. In the result, the es- 
timate for the critical length scale L cr it gives moderately 
higher values L cr it oc (l + z s ) 01 , where z B denotes the time 
of shock formation. This means that if shocks happen at 
higher redshifts, also the upper bound for the mass of the 
halos which can still be provided with dense cold gas via 
the attached filaments remains almost unchanged. We con- 
clude, that for halo masses exceeding the above given range 
3 x lO n Af < A/haio < 10 13 M Q (1 + z s ) - 3 a supply of the 
inner halo with cold gas will be terminated. 

(v) Observations: The temperature of the shocked gas 
in the filaments scales as oc L 2 . Thus, for perturbation scales 
of L w 10 Mpc/h the temperature of the filament gas reaches 
almost 10 K. Those filaments can be attributed to the bulk 
of the WHIM as predicted by hydrodynamical cosmological 
simulations. The observation of WHIM filaments at these 
scales remains challenging. The core region of the filament 
becomes too small and the probability that a line of sight in- 
tersects with this particular region becomes extremely small, 
as well. Thus, the filaments are hot almost everywhere and 
do not exhibit significant variations which could serve as a 
characteristic tracer for observations. However, if by chance 
the line of sight hits the core region then absorption by hy- 
drogen could be observed at high column densities. Since the 
spatial extension of the core scales as ex L -3 ' 57 and the den- 
sity as oc L 5 ' 7 , the column density of neutral hydrogen scales 
as oc L 7,83 as long the core is still optically thin. Using the 
values for L = 4 Mpc we obtain an estimate for the column 
density AT HI = n H i A D « 2 x 10 14 (L/4Mpc) 7 ' 83 (1 + z s f . 
For 2^3 and for L ^ L CT n the assumption of an opti- 
cal thin core holds. At higher redshifts or for considerably 
larger L the core may become shielded against the UV radia- 
tion. In this case, radiative cooling is not longer balanced by 
UV heating and further contraction will occur. However, at 
some stage the medium becomes optically thick with respect 
to recombination radiation. The radiation is not longer able 
to leave the cold core freely, but interacts significantly with 
the gas before leaving the core region. Thus, for a primordial 
matter com position, the cooling becomes less efficient (see 
iGnat &: Sternbe rg 2007) and for a full description radiative 
transfer effects must be taken into account. In this context, 
and as a a goal for future studies, it is needed to consider 
possible ionization states and the corresponding column den- 
sities for heavy elements. Temperature and density profiles 
with partly steep gradients may cause characteristic absorp- 
tion features. The main uncertainty, however, remains the 
abundance of the heavy elements in the filaments far from 
halos/galaxies. 
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